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ABSTRACT 


Experimental data are presented for the growth of thrombi 
(blood clots) in a stagnation point flow of fresh blood. Thrombus shape, 
size and structure are shown to depend on local flow conditions. The 
evolution of a thrombus is described in terms of a physical model that 
includes platelet diffusion, a platelet aggregation mechanism, and 
diffusion and convection of the chemical species responsible for 
aggregation. Diffusion-controlled and convection-controlled regimes 
are defined by flow parameters and thrombus location, and the 
characteristic growth pattern in each regime is explained. Quantitative 
comparisons with an approximate theoretical model are presented, and 
a more general model is formulated. 
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I. INTRODUCTION 


The coagulation of blood is a complex phenomenon involving bio- 
chemistry, surface chemistry, and fluid mechanics. The purpose of this 
experimental and theoretical research is to elucidate the influence of the 
local flow conditions of the blood on the formation of a thrombus, particu- 
larly on artificial surfaces. To give an example of the role played by the 
flow, first consider the case of stagnant blood in contact with a foreign 
surface such as a test tube. The blood solidifies into a red clot by the 
formation of a stringy fibrin mesh which entraps the other constituents of 
the blood. In contrast to this, when flowing blood contacts a foreign surface 
or injured tissue, white solidified masses are observed. These masses, 
called thrombi, consist largely of platelets which are cells roughly one 
micron in size and occupy about 0. 1% of the volume of whole blood. The 
presence of flow indeed has a marked effect on the mechanism by which 
blood solidifies in response to external stimuli. 

It is generally accepted that high flow rates and high values of shear 
at the surface (short of that necessary to damage red cells) tend to minimize 
the incidence of thrombosis. Large shear forces on blood elements loosely 
attached to the surface may detach these elements before an actual throm- 
bus is formed. A second effect is that high flow velocities imply shorter 
times in which the blood is in contact with the foreign surface. Some chemi- 
cal reactions are required either to cause platelet aggregation or to activate 
the coagulation sequence. The shorter the residence time of the blood near 
the surface the less time is available for these reactions to proceed and, 
therefore, one may expect that less thrombus or no thrombus will form 
under conditions where the residence time is short. On the other hand, a 
mechanism exists by which increased flow tends to lead towards greater 
thrombus buildup. In the presence of high surface shears the rate at which 
platelets diffuse to the surface is increased. Therefore, if the surface is 
such that platelets stick on arrival at the surface, a high flow rate would 



lead towards a more rapid accumulation of platelets and a rapid growth of 
the thrombus. Thus, increased flow could lead to either increased or 
decreased thrombosis, depending upon which of the above mechanisms is 
dominant. As will be described later, some evidence for each of these 
mechanisms has been observed. 

Over the last several years a research program has been carried 
ouF ; to study these fluid mechanical effects both experimentally and 
theoretically. A carefully controlled stagnation point flow experiment has 
allowed the direct observation of the above factors involved in the growth 
of thrombi on a foreign surface. The experimental observations suggest 
that diffusion of platelets to the surface, attachment of white cells within a 
well-defined radius of the stagnation point, and subsequent platelet aggrega- 
tion in the same region are important flow dependent processes at low flow 
rates. However, at sufficiently high flow rates platelet aggregation and 
thrombus growth occur only at surface imperfections and tend to assume a 
characteristic wedge shape with time. 

A series of observations is presented of the growth history of indivi- 
dual thrombi at various radial locations and flow rates. The shape, approxi- 
mate thickness, and internal structure have been monitored over a sufficient 
range of flow conditions to demonstrate the variety of possible flow dependent 
growth patterns. A physical model and scaling parameters for the spatial, 
and temporal behavior of a thrombus is also presented. A "complete” theo- 
retical model, the final objective, remains to be developed. However, sev- 
eral clearly important features of the growth process have been modeled 
and approximate solutions obtained. Finally, quantitative comparisons with 
experimental data have been made where possible, and the more general 
analytical model outlined. 

The evolution of a thrombus on an artificial surface is discussed in 
Section H. The experimental equipment and procedure is described in 
Section III. The theory of a stagnation point non-Newtonian boundary layer 
is reviewed in Section IV. In Section V the detailed experimental data are 
presented for the change of shape, size and structure of thrombi as a func- 
tion of time. A general formulation of the theoretical problem is given in 
Section VI and the scaling laws derived from these equations are compared 
to the data. Section VII is the concluding discussion. 
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H. THE EVOLUTION OF A THROMBUS ON A FOREIGN SURFACE 


When fresh blood is brought into contact with an artificial surface, 
a very complex chain of events is triggered which differs in numerous 
respects from the interaction of blood with its natural container, the blood 
vessel waits. In the absence of any disease which might predispose the 
blood to some type of hypercoagulability, the walls of the blood vessels have 
an antithrombogenic influence on the blood that designers of artificial devices 
would be pleased to emulate. When this cell wall (the endothelium) is injured, 
physically or chemically, a chain of events occurs which may lead to the 
growth of a thrombus. Although there are similarities to growth on an 
artificial surface, the differences, which have to do with agents released 
from the injured tissue, are sufficiently profound that the reader is cautioned 
not to directly extrapolate our description of the observed growth process to 
that occurring on an injured cell wall. 

The first microscopically observable occurrence in the sequence is 
the adhesion of individual platelets to the surface 1 . This does not occur 
immediately, however, and the small delay time of the order of a minute 
is thought to be the time required for a layer of protein to be adsorbed onto 
the surface from the blood plasma and modified so that it becomes attrac- 
tive to platelets. Although the type of surface must have some effect on the 
manner in which this layer is deposited, the exact nature of the interaction 
is unknown and a matter of controversy. A recent study by Friedman, 

(7 ) 

et al concludes that the platelet deposition rate is not a strong function of 
the type of surface, which we also noticed in earlier experiments. This 
could be the result of the adsorbed protein layer which physically separates 
the surface from the blood. 

The second step is the deposition of a monolayer of platelets on the 
surface. The method of adhesion of the platelets to the protein layer is not 
known. When the process stops at this step in the sequence, which we have 
observed often under some flow conditions, no thrombus forms during the 
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duration of the experiment 

The third step is platelet aggregation. During this phase the platelets . 
already on the surface are activated in some manner which causes them to 
move together along the surface into clumps. This may be accompanied by 
distortion of the platelet cell membrane, e. g. , the platelets swell and extend 
outward projections called pseudopodia^ 8 ' The peculiar behavior of plate- 
lets during aggregation has been studied a great deal. See, for example, 

the other. articles in the collections edited by Sherry, et al^ 10 ’ ^ and by 

( 12 ) 

Johnson and Guest . Other books of interest on the subject are those by 

Marcus and ZuckeK * and the books edited by Kowalski and Niewiarowski^ 14 ^ 

(15) 

and by Johnson, et al 

These platelet aggregates are illustrated in a series of scanning 
electron micrographs in Figs. 1-3. The first shows a clear dividing line 
between a region of heavy platelet aggregation (aggregates in white) and a 
region which is comparatively unaffected. The second micrograph. Fig. 2, 
shows individual platelets which have sent out tendril- like pseudopodia but 
have not yet aggregated. The larger ( ~ 8 /x ) smoother objects are red cells 
which were not actually attached to the surface at this point in the experi- 
ment and the presence of which is misleading. They appear only because of 
incomplete washing of the coverslip. In Fig. 3, the platelets are no longer 
individually recognizable and have aggregated into much larger clumps, the 
surfaces of which appear almost hairy. 

Since we do not see platelet aggregates everywhere on the foreign 
surface, one might say that the aggregation cannot be simply a matter of 
reactivity between the surface and the platelets. However, we observe that 
surface imperfections or inhomogeneities (which may be seen optically) 
often are the sites of significant platelet aggregation. Complementary to 
this is the observation that physical bumps in the surface, but of the same 
general material, may significantly disturb the flow and yet not cause aggre- 
gation of platelets. One is forced to conclude that some property of the 
surface imperfection, whether it be a change in electronegativity or surface 
energy or whatever, acts as a triggering mechanism. 

The biochemical processes which lead to platelet aggregation have 
been studied extensively. Most of these studies are performed on platelets 
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Scanning electron micrograph of an aggregated and unaggregated 
region. 




Fig. 2 Scanning electron micrograph of the unaggregated region of 

Fig. 1. Note the density of the individual spread platelets with 
minimal interaction between them. 
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Scanning electron micrograph of the aggregated region of Fig. 1. 
Note the density of platelets and small platelet aggregates 
between the larger major aggregates. 


Fig. 3 


suspended in plasma. Numerous chemicals which lead to aggregation have 
been identified^ \ These studies do not, however, define which chemical 
species actually play the dominant role in platelet aggregation in flowing 
blood. One of the most likely candidates is ADP (adenosine diphosphate). 
Small concentrations of ADP lead to reversible platelet aggregation, i. e. 
the platelets are observed to aggregate into clumps and then deaggregate into 
individual platelets. Larger concentrations of ADP lead to apparently per- 
manent aggregates. This process is believed to be accompanied by release 
of further ADP from the platelets themselves. The rate of accumulation of 
platelets on the surface is highly dependent on the local velocity profile 
which, as will be shown later, determines the relative diffusion rates of 
platelets and the species responsible for aggregation. 

The processes of thrombus formation and blood coagulation (forma- 
tion of a fibrin and red cell mass), while distinct, may not be entirely inde- 
pendent. Thrombin, which is one of the precursors in the chain of reactions 
which lead to fibrin formation, is also known to lead to platelet aggregation. 
Furthermore, the later stages of thrombus formation involve the develop- 
ment of a fibrin mesh in the stagnant regions between platelet aggregates, 
and the resultant trapping of red cells within the fibrin (see Fig. 4). 

For the purpose of the present discussion, the precise nature of the 
chemical species which leads to platelet aggregation is not important. We 
shall assume that the surface imperfection which initiates the thrombus 
causes local platelet aggregation. Our primary concern will be with the 
influence of flow on the further development of the thrombus. In attempting 
to analyze this we will further assume that the aggregated platelets release 
some species which can cause further aggregation as it is convected down- 
stream and diffuses away from the streamlines. While future work may 
reveal the exact chemistry of the process, it is possible to explain some of 
the observed features of thrombus growth based on the fluid mechanics, 
requiring only the above simplified assumptions concerning the complex 
biochemical phenomena. 
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Fig. 4 Scanning electron micrograph of thrombus in an advanced stage . 
Note the red cells trapped within the fibrin strands. 
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III. THE STAGNATION POINT FLOW EXPERIMENT 


In this section, the experimental procedure used for the thrombus 
formation study will be described. The general procedure is now routine 
and is utilized to study thrombus formation both as a function of hemo- 
dynamics and material surface properties. 

Mongrel dogs weighing approximately 20 to 25 kg are anesthesized 
by intravenous injection of sodium pentobarbital. The right cartoid artery 
of the dog is exposed and everted over a tube which is then connected to the 
flow chamber. Seven flow chambers have been designed for the experiment 
which have inviscid flow parameters varying from 3. 8 sec’^ to 85 sec" 
Figure 5 is a schematic diagram of a flow chamber. The chamber is 
designed such that the first foreign surface which the blood contacts is the 
surface being studied. For this series of experiments, polyurethane 
surfaces were utilized. Polyurethane was selected because of its optical 
and relatively non-thrombogenic properties. It was solvent cast on 25 mm 
diameter microscope coverslips and baked in a low temperature oven with 
nitrogen gas flowing over the surface. This' baking procedure firmly bonds 
the polyurethane to the glass surface, thus preventing detachment during 
the experiment. 

The chamber is initially primed with saline to prevent the formation 
of an air-blood interface and the trapping of air bubbles in the chamber 
which result in flow pulsations. The blood is allowed to flow through the 
chamber at a flow rate of 2 cc/min and is collected in a withdrawal chamber. 
The blood is not returned to the animal, thus preventing any contamination 
of the blood contacting the surface by any previously "activated" blood. 

The flow rate through the chamber is controlled by a withdrawal system. 

The system is designed such that pulsations are minimized. For this 
series of experiments the animal was not heparinized. Heparin is added 
to the blood as it flows from the flow chamber into the withdrawal system. 
This prevents clotting in the withdrawal lines and collecting chamber. 
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Fig. 5 


Schematic diagram of flow chamber in the stagnation point blood 
flow experiment. 


During the experiment, the blood-polyurethane surface is continually 
observed by dark field reflected microscopy. The location of the stagnation 
point and any thrombi which form are noted. When a thrombus is observed, 
its growth is documented by visual observations, 35 mm photographs, and 
16 mm movies. Time at which flow distortion begins to occur, reverse or 
lateral flow, and red cell trapping are observed visually.. Approximate 
a gg r egate depth is obtained by utilizing a 22 X objective lens and a micro- 
meter scale on the focus adjustment of the microscope. The lens is first 
focused on the blood-polyurethane surface, and then focused through the 
aggregates on the blood cells flowing over the aggregates. 

Local white cell velocity data and thrombi growth profiles are 
obtained from the 16 mm movies. Profiles and interaggregate distance are 
also obtained from 35 mm slides. 

At the conclusion of the experiment, the dog is sacrificed and the 
coverslips are stained with Wrights stain and examined. 
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IV. THE STAGNATION POINT NON-NEWTONIAN BOUNDARY LAYER 


The flow of blood on the surface of the coverslip is modeled as 
axisymmetric stagnation point flow of a continuous viscous fluid having a 
non-Newtonian stress-strain relationship. In the absence of viscous 
effects, the velocity components would be given by the potential flow solution 

u = +Ar, v = -2Az (1) 

where r, z are the radial and normal coordinates and u, v the corresponding 
velocity components. The flow field, is illustrated schematically in Fig. 6. 

The constant A is the stagnation point flow parameter, which represents the 
radial velocity gradient of the fluid at the edge of the boundary layer as it 
accelerates away from the stagnation point. The above velocity field is only 
valid for distances from the surface large compared to the boundary layer 

/ I 

thickness. For a Newtonian fluid this thickness is given by'’ 

6 n = 2.4(v/A) l/2 (2) 

2 

where vs. 04 cm /sec is the kinematic viscosity of blood at high rates of 
shear where blood is Newtonian. To get some appreciation of relative sizes 
we can calculate that for A = 10 sec this thickness is 6^ = 0. 15 cm which 
is about the characteristic size of our flow chambers and is about 190 red 
cell diameters. 

It is in fact the presence of the red cells which causes blood to have 
a non-Newtonian behavior at low shear rates, for the disk-like red cells 
have a propensity to form long chains (rouleaux) not unlike stacks of poker 
chips. The adhesion of the red cells to one another, which is thought to be 
caused by the adsorption of fibrinogen on the surface of the cells, may be 
broken at higher shear rates and therefore blood is what is termed in 
rheology a pseudoplastic, or shear-thinning fluid. A non-Newtonian stress- 
strain law which adequately fits the experimental data v 1 is the power law 
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Fig. 6 
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Schematic diagram of a stagnation point flow field, showing 
non-Newtonian velocity boundary layer and platelet diffusion 
layer. 
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(3) 


T = K( 9u/a z) N 

where r is the shear stress, 9u/9z the rate of strain and K, N are 
constants. Blood of normal hemocrit (45% red cells by volume) may be 
represented^ 1 ^ by the values N = 4/5, K = 0. 1 in cgs units. Cor- 
respondingly, one may consider blood to have a rate of strain dependent 
viscosity of p, & 0. 1 (9u/9z) 2 gm cm 1 sec 1 . The above approxima- 

tions should be adequate for the boundary layer analysis which follows. 

The axisymmetric stagnation point boundary layer equations are 

pu 9 u/9 r + pv d u/d z = dr / d z ~ dp/d r (4) 

9 u/9 r+u/r+ 9 v/9z = 0 (5) 


for the conservation of radial momentum and conservation of mass. The 
shear stress is given by Eq. (3) and the pressure gradient is 9p/9r = -p A r 
from the irrotational flow represented by Eq. (1). It is convenient to non- 
dimensionalize all the lengths and velocities by the reference length 
£. - (v -gj A) 1 / 2 and the velocity U = (v -^A) 1 / 2 where is a reference 
viscous diffusivity to be chosen momentarily and A is the stagnation point 
flow parameter. Letting 


(u, v) = (u/U, V/U); (r, z) = (r/i , z/l) 
we have for the dimensionless conservation equations 

N- 1 


u 9 u , - 9u - 
+ v 


9 r 


9 z 


■(H) 


9 2 u 
9 z 2 


9 u 
9 r 


+ — + 


9 v 
9 z 


= 0 


if the reference kinematic viscosity is chosen as 

N-l 


v 


NK A 


R 


+ r 


( 6 ) 

(7) 

( 8 ) 

(9) 


- 17 - 



The proper nondimensionalizing length and velocity are then 


fNK A 


N-2 


\ 1/2 / 

) 


’) 


1/2 


P / \ P 

The above system of partial differential equations, (7) and (8), may be 

reduced to one ordinary differential equation in one independent variable 

using the similarity transformation implied by the following equations 

1-N 


( 10 ) 


q = z r 


1+N 


(II) 


u = r $ 1 (r) ) 


( 12 ) 


-z 


3N+_l\ 4_ / 

N + 1 J 7] + ^ 


1-N 
1 + N 


r 


The differential equation is 


(13) 


, ,N-1 , , . , 3N + 1 , , . ? 

^ +HnTTM -W ) + 1 = 0 


(14) 


with the boundary conditions 


<P (0) = (£ * (0) = 0, <f>' («) = I 


(15) 


The numerical solutions to (14) may be found in Ref. 20 where <£'(q) and 
<j>' ' (0) are tabulated fox various values of N. 

The velocity gradient and the shear stress at the surface are given 
in physical coordinats for a non-Newtonian fluid by 


1 


S u\ / pA 3 r 2 \ 

d z ) \ NK / 


(0) 

(16) 

4 X ’ \ 




/ . 3 2 \N+1 




T - K j P A r ] 

w l NK J 

to' 

o' 

(17) 


The relevant constants that have not yet been given are T (0) « 1. 25 for 
any N near unity and p « 1.06 gm/cm . 
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For a Newtonian fluid, the shear at the surface can be shown to 
increase as the 3/2 power of the flow parameter, A, and to increase 
linearly with distance from, the stagnation point. Thus, for a given 'flow 
parameter, the shear is zero at the stagnation point and increases as one 
moves outward. The shear at any given distance from the stagnation 
point is larger the higher the flow parameter. The radius at which a given 
shear is reached decreases as the 3/2 power of the flow parameter. Thus, 
for a situation in which local shear is the significant parameter, one would 
expect similar phenomena to occur at radii which decrease with increasing 
flow as the 3/2 power of the flow parameter. For some of the results to be 
described later, the departure of blood from a Newtonian viscosity Is 
significant. In this case the shear still increases with radial distance from 
the stagnation point. However, its variation with distance is different. It 
can be shown, however, based on an analysis of non- Newtonian stagnation 
flow, that the radius at which a given shear is reached still scales as the 
inverse 3/2 power of the flow parameter. 

For comparison to experiment it is useful to approximate the velocity 
profile in this thin viscous layer as a linear variation with the distance z 
away from the surface. In the viscous layer then 

u = p rz (18) 

where we have defined a convenient parameter p = (3u/dz) r - . 

The parameter |3 is theoretically well defined only for the ideal 
stagnation point flow, the streamlines of which are hyperbolic in shape as 
shown in Fig. 6. Since the flow in the experimental flow chambers does 
not correspond exactly to this ideal, the parameter P has been measured 
empirically. This also allows the determination of an effective inviscid 
stagnation point flow parameter, A, for a given flow rate of blood through 
a given chamber. 

This discussion of the velocity field in the flow experiment forms a 
framework for the following presentation of the data taken on the growth of 
thrombi. 
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V. THE GROWTH OF A THROMBUS: EXPERIMENT 


Numerous experiments have been carried out according to the 
procedure outlined in Section III. The data to be presented here have been 
selected as typical of the numerous thrombi we have observed. Other data 
were rejected for various reasons, e. g. , some thrombi could not be 
observed over a long enough time interval or over a large enough area; 
others occurred so close together in one run that their growth was dom- 
inated by mutual interference effects. The four thrombi (Labeled A, B, C, 

D) which we show here are the best documented examples. 

The photographs of the four thrombi are black and white reproductions 
of single frames of the color movies taken through the optical microscope. 
The first photograph. Fig. 7, shows thrombus B at a very early stage. 

What can be seen as a faint airfoil shape is the region of platelet aggregation 
and white cells in the wake of the point of initiation. A schematic picture of 
the growth of the region of platelet aggregation with time is furnished in 
Fig. 8. A picture of the same thrombus at a later intermediate time is given 
in Fig. 9. Note the swept back wedge-like shape of the thrombus and the 
mottled appearance and nonuniform nature of the platelet aggregation. The 
picture of thrombus A in Fig. 10 clearly shows the characteristic shape of 
the head or leading edge of the wedge thrombus, and the initiation point 
(a surface imperfection) shows up as a bright spot. The last photo of throm- 
bus A in Fig. 11 shows the shape at a later time when the growth has been 
greatly influenced by the radial nature of the flow. The radius of curvature 
of the indentation of the head of the thrombus is nearly equal to the radius 
from the stagnation point. Some thrombi, such as thrombus D, eventually 
grow to envelop the stagnation point, and thereafter rapidly cover the entire 
surfac e. 

Quantitative data chronicling the growth history of the thrombi are 
presented in three ways. First, in Figs. 12-15 the outlines of the greatest 
extent of the platelet aggregation as perceived visually are given as a 
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Fig. 7 Thrombus "B" at 21 min. 00 sec. after start of the experiment. 
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Fig. 8 Schematic diagram of thrombus growth. 
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Fig. 9 Thrombus "B" at 36 min. 30 sec. after start of experiment. Note 
this photograph is a composite of two photographs taken during the 
experiment, in order to show the increase in length and width of 
the thrombus. 
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Fig. 10 


Thrombus "A" at 1 9 min. 


00 sec. 
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Fig. 11 Thrombus "A" at 38 min. 30 sec. 


26 













function of time for each of the four thrombi. Secondly, data on the growth of 
the head toward the stagnation point and on the width of the head is presented 
in Figs. 16-19. These measurements are illustrated in the schematic of 
Fig. 8 . In Fig. 20 we have graphed the thrombus head position as a function 
of time for the four thrombi grown under different flow conditions. The 
third set of data is given in Fig. 21 and represents statistics on the inter- 
aggregate spacing in a thrombus as a function of distance from the stagnation 
point and time. The data is presented in the form of histograms of measured 
interaggregate distances at three times during the experiment. A "best fit" 
curve is drawn through the total set of measurements at each time. 

It is clear that a theoretical description of these effects must con- 
sider the transport to the surface of molecular and cellular matter suspended 
in the blood as well as the surface reactions. The characteristic size, shape 
and structure of the wedge thrombi discussed above are qualitatively well 
understood in terms of the model which is formulated in the following section. 



Ar $ W/2,mm 



Fig. 16 Upstream growth of the thrombus head, A r, from the thrombus 
initiation point, r^, and the half width, w/2, of the thrombus as 
a function of time from the initial contact of the blood with the 
surface. Thrombus "A”; r^ = 1.40 mm, A = 16.4 sec - -*-. 
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Fig. 20 Growth of heads of four thrombi with time. 
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Fig. 21 


Total number of aggregate pairs vs. interaggregate distance 
A = 10.4 sec" 1 , rj = 0.9 mm. 
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VI. THE GROWTH OF A THROMBUS: THEORY 


Part I - Early and Intermediate Time 

In the section describing the evolution of a thrombus, a two-part 
mechanism of platelet aggregation and accumulation is postulated which 
utilizes the concept of an activating species which is contained in the 
platelets themselves. When a platelet encounters a critical concentration 
of this species near the surface, it sticks to the surface and joins in the 
aggregation, in which process more of the species is released to activate 
other platelets. This process is graphically described in the flow chart of 
Fig. 22. The growth of the aggregate can be limited only by a shortage in 
the release rate of the activating species or a shortage in the arrival rate 
of platelets. Now, we would like to quantify these concepts in the con- 
struction of a mathematical model capable of describing the features of 
thrombus growth which have been documented in Sec. V. 

The platelets are the dominant constituent in the thrombi we observe 
We shall derive a conservation equation for them presently. Before doing 
this, however, we must give attention to the dominant transportation mech- 
anism for such macroscopic particles. In Ref. 3, several mechanisms for 
platelet transport to the surface are considered. It is shown that direct 
impingement of the platelets as they are carried by the flow to the surface 
along streamlines is a negligible contribution to the platelet arrival rate 
and is inadequate to explain observed platelet growth. Similarly,- because 
the platelet is so large compared to molecular species, the diffusion by 

Brownian motion is very small. The random diffusion coefficient for plate- 

~92 ~52 

lets is roughly 10 cm /sec. compared to ~10 _ cm /sec for low molecula 
weight species, and hence the flux of particles by ordinary diffusion is very 
small. 

The important transport mechanism -is redcell augmented diffusion. 
The shearing motion of the velocity gradient in the boundary layer causes 
the red cells, both as individuals and in rouleaux, to tumble at a rate pro- 
portional to the velocity gradient. The velocity perturbations caused by 
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RATE OF NEW PLATELETS 


Fig. 22 Platelet aggregation sequence. 
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these rotating disk-like cells result in a significant enhancement in the 
mixing rate, not unlike turbulence in nature. The essential effect of the 
velocity perturbations induced by the red cells rotating in response to the 
velocity gradient is to make the augmented diffusion coefficient also pro- 
portional to the velocity gradient. 

The first order effects of red cell tumbling will be estimated with 

a fairly elementary model. Due to the velocity gradient '( 9u/8 z) in the 

boundary layer, a red cell rotates about its diameter with an average 

1 0 u 

angular velocity CO = , as required by a zero net torque' on the 

particle. The velocity field consists of two components; one due to un- 
bounded shear flow, and one due to particle rotation. Representing the 
rotating red cell as a sphere of radius b, the latter perturbation decreases 
with distance^ 1 ^ as ~ (b/r)^, whereas the former varies as ~(b/r)^. For 
dilute suspensions (low hematocrit), the only long-range effect is due to 
red cell rotation, and the effects of red cell interactions are also negligible. 
Completely random interactions between the small particles (platelets) and 
the red cell velocity fields will be assumed, with a characteristic frequency 
f and "mean free time" t given by the relative velocity Au and average dis- 
tance 2Tr between the red cell and particles; 


f = — = ^u 
t Ar 

hence t = ■ = — 1 — 

Zc d Az 2 co 

. J 

Furthermore, A^ ~ l/2 C ' , where C is the number density of red cells. 

Finally, the characteristic velocity perturbation is q'~cob (b/ZAr) . Now, 
for a "random walk" process, the "mean free path" A. — (Dt)^^ where D is 
the diffusivity. Since X~q’t as well, 

2 

D~(q , )t~8W 2 b 2 (b/2A'r) 4 
or D ~4b 6 C 4 / 3 8u/3z 
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Finally, in terms of the volumetric concentration C , 

D ~ b 2 C^ 4 / 3 9u/8z (18) 

■with the constant of proportionality expected to be of the order one. 

The above result is an estimate which is valid, at best, at low 

hematocrit, for unbounded flows, and small red cell Reynolds number. 

In Ref. 21 it is shown that this derived diffusion coefficient is a 

good approximation to data taken in a simulation experiment using neutrally 

buoyant spheres in a pipe flow geometry even at moderately high particle 

concentrations, velocity gradients and particle Reynolds number. 

At very low shears this type of diffusion coefficient is inappropriate 

because the red cells form rouleaux, i. e. , they cling together in long chains 

resembling a stack of poker chips. The diffusion coefficient varies as the 

square of the largest radius of the particle, so that an increase in particle 

(rouleaux) length with decreasing shear will reduce the dependence of the 

diffusion coefficient on the shear. 

{ 22 ) 

A model v ' of rouleaux formation that appears to correlate blood 
viscosity data at low shear yields the following expression for the length to 
diameter ratio J: 


J « 1 + 


M) 


- 1/2 


(19) 


Thus, as 9u/9z— ►O, J — (9u/9z) ^/ 2 and the proper length in Eq. (18) is 

(a) 


J b' 


A diffusion coefficient limit at very small shear rates is then obtained: 


D — - 2 C 4 / 3 ( Jb) 2 ~3 C 4 / 3 b 2 
P 8z p 

For normal blood, then 

D -*1. 5 x 10 ^ as 3u/9z — 0 
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As shear rate increases, J 1, and the previous result is obtained. A 
reasonable matching point is at 8 u/9 z = 1. 0, so that 


and 


D « B (9 u/9 z) M 
B = 1. 5 x 10' 7 


M = 


I 


/ 


0 for 9 u/9 z <1 


1 for 9 u/9 z > 1 


( 20 ) 


The equation expressing the conservation of platelets is the usual 
steady state boundary layer species conservation equation, which is written 
as 


9 C 

u{r, z) 8 J + v(r, z) 



9 

8 z 


( 


■ 9 C 



( 21 ) 


where C p is the volumetric concentration of platelets and u, v are the local 
velocity components in the momentum boundary layer which are given in 
Sec. IV. Azimuthal and radial diffusion terms have been dropped as being 
small compared to the normal diffusion, and the term involving an azimuthal 
velocity will not be present until the thrombus is large enough to appre- 
ciably distort the flow. There are no production or loss terms in this equa- 
tion because the only loss of platelets is at the bounding surface. The outer 
boundary conditions are that the concentration of platelets assumes the free 
stream value on the axis and at the edge of the boundary layer 


C p (r=0, z) 



lim 

zlarge 


C 

P 


(r, z) 



( 22 ) 


The surface boundary condition we shall leave for later consideration. 

In Ref. 2, a similarity solution to Eq. (21) is obtained in order to 
get the flux of platelets driven to the surface by the augmented diffusion. 
This flux can be shown to be, for a Newtonian boundary layer and M = 1. 
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F = 
w 


D 

P 



platelets 

2 

cm sec 


= 2x10 


-4 


(A sec 


5/3 


(r, cm) 


2/3 


( 23 ) 


The rate of growth of the thrombus in the z direction can be estimated 
dli 

then to be ■■ = (flux of platelets to thrombus) (volume of platelet) 


= F (4x10 12 cm 3 ) 
w ' 


(24) 


-(8x10 A 3 /” 3 (r, cm) 2 / 3 C cm/sec 


For thrombus C, and A= 17 sec -1 , r = 0. 155 c, and C = 4. 6 x 
10 cm . The rate of increase of height of the thrombus is 5 x^O - ^ cm/sec 
or 3 microns/min. The estimates of observed thickness indicate a smaller 
rate of growth, like one micron per minute. The halfwidth of thrombus C 
grew at a rate of about 15 microns/min, however. (See Fig. 18). Since 
the growth rate calculated should be the maximum rate of accumulation of 
platelets, it is apparent that the rate of platelet aggregation is probably not 
limited by platelet arrival rate but rather by the diffusion of activating species. 

The equation for the conservation of an activating species has a simi- 
lar form, only now the diffusion in each of the three directions can be of the 
same magnitude, so the equation has the form. 


3 C 8 C 9 C 

0l + u(r,z) -g / +v(r,z) - a 



1 _ 9 

r - 9 r 



2 

9 C 


9 z 


( 25 ) 
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where C is the volumetric concentration of activating species and D is the 

** ct 

corresponding diffusion coefficient, which will be taken to be a constant. We 
may not, in general, omit the time dependent term here because the concen- 
tration contours follow the thrombus surface quite closely and the thrombus 
growth time history is what we seek. The outer boundary conditions for C 
are that vanishes for large distances from the thrombus in each direction. 
The initial condition is that vanishes everywhere initially, except at the 
initiation point of the thrombus (r = r., 6 = 0), at which point it has some 
given value which is greater than a critical concentration for platelet aggre- 
gation. 

The phenomenon of the thrombus exuding a particular species which 
in turn diffuses into and is convected by the boundary layer flow has several 
interesting features. One can best discuss the various regimes of interest 
in a coherent manner by reference to the Thrombus Growth Map in Fig,, 23. 
The various shapes of thrombi in a stagnation point flow are illustrated. 

For early time and for radii of the same size or smaller than a diffusion 
length, the thrombus is diffusion dominated and will tend to have near equal 
lateral and radial dimensions. For large flow parameters or for points of 
initiation large compared to a diffusion length, the thrombus is convection 
dominated and appears to have a very swept-back shape. 

We can get an estimate of the velocity of the thrombus edge growth 
(the moving platelet aggregation front) in the following way. -Assuming that 
the thrombus is exuding the activating species, at some time r the locus of 
the contour of constant critical concentration of this species is assumed to be 
at a distance 5 from the current thrombus position. The thickness of this 
diffusion layer must depend on a balance between the diffusion perpendicular 
to the layer and the velocity component perpendicular to the thrombus edge. 

A good estimate of the size of this diffusion layer is 


6 


D 

a 

u(r, z) cos a 


( 26 ) 


where a is the angle between the direction of the boundary layer velocity 
vector and the thrombus edge normal. There is a question as to what 
vecocity to use here, i. e. , how far up into the velocity boundary layer does 
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Fig. 23 Wedge thrombus growth map. 
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this diffusion process extend? It must be of the same order as the sideways 
diffusion, so we evaluate the velocity in the boundary layer at a height equal 
to 6. We shall use 


u(z = 6) 


8 u _ Ar 4>" (0) 5 

8z "vwx- 


(27) 


for the velocity in the boundary layer in the previous formula for diffusion 
thickness. Solving for 5 then gives 


6 = 



D 

a 


cos a 



D ^ 
a 


1/2 


1/2 


11 (0)r cos a 


(28) 


This is a measure of the distance in which the diffusion front moves 
/ a u \ - 1 

in the time I - q- — — cos a J and defines the region within which the concen- 
tration of activating species is above some critical level. This region will 
march forward in time as platelets arrive within the region, stick, and 
release more of the species. We can estimate the characteristic time r 
for the surface to collect a particular surface density of platelets by making 
the ratio of the flux to the wall of platelets to the surface density of platelets. 
This surface density should properly be that necessary to release the criti- 
cal concentration of the activating species. This would be a very interesting 
number to derive from the measured growth velocity data once this simple 
model is complete. 

Then 


T 


F 

w 


N 

pw 


d a c 

_£ £ 

N 9 z 

pw 


-1 

sec 


( 29 ) 


is a fair estimate of the characteristic time to deposit a particular surface 
density of platelets. 
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The velocity of the edge of the thrombus is the distance in which the 
region 6 of sufficient activating species advances as it fills up with platelets 
in time r . This speed is 


u 


th 



D 

P 

N 

pw 


8 C / D 

EL / a 

9 z i 9 u 

\Tz cos ° 



(30) 


In the Newtonian region where the platelet diffusion is augmented M = 1 and 
the flux is 


F = D 
w p 


3 C 

P - 





(31) 


So the thrombus front velocity is 


A 3 / 4 

U th “ 4 


N 


^ (bV/ 2 ) 


1/3 


pw 



D 

a 

<t > 11 ( 0 ) cos 


a 



(32) 


The dependence of the growth velocity- on the parameters is very interesting. 
The strongest dependence is on the flow parameter and on the platelet free 
stream concentration with lesser dependence on. the species diffusivities and 
radius. 


The required surface density of platelets to get a critical concentra- 

C * ? 

tion of activating species is taken as N = 1. 5 x 10 platelets/cm . Fried- 
(7) P w 5 t 

man v , for example, quotes about 70 x 10 plate lets/cm total surface 

density deposited after six minutes or so. These numbers are consistent 

in the sense that the number of platelets necessary per unit area to continue 

the growth of the aggregation region is likely to be smaller than the total 

surface density already there. 

Note that the geometrical factors in the growth velocity, u , ~ r^^* 

~l/Z . . ' 

cos a , would indicate fastest growth where the edge is aligned with the 

flow ( large a ) and at larger radii, which is consistent with our observations. 
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The expression for the thrombus growth velocity' in (32) should 
collapse the data taken for the four thrombi at different" flow conditions. 
Taking u ^ = dr/dt and integrating to find the front position versus time, 
r(t), we get 


(r 


5/6 5/6, 1/2 

' - r. ' ) cos ' 

l ty 


5 

24 


C B 2 / 3 
. P=o 


N 

pw 


Da 

* 1/2 * 


1/2 

~ ) (33) 

( 0 ) 


The factors on the right hand side of the equation are fixed constants for any 
blood sample. The factors on the left side are the scaling factors for throm- 
bus position as a function of time and angular position for different flow 
parameters. 

In Fig. 20 we have graphed the thrombus head position as a function 
of time for the four thrombi grown under different flow conditions. The 
agreement is excellent considering the uncertainties in the concentrations 
and diffusivities . 

Part II: Late Time 

The asymptotic structure of a growing wedge thrombus is characteri- 
zed by a large number of well define platelet aggregates. As seen in Figs. 

9, 10 and 11, the aggregates form in the w.ake of the initial disturbance 
(thrombus head) and grow with time. As they grow, their total number 
decreases, and the interaggregate distance increases. Thus, the ratio of 
aggregate size to separation distance appears to be fairly constant with time, 
and the characteristics aggregate dimensions apparently "scale" with throm- 
bus size. Since aggregate growth has been shown to be due primarily to 
platelet diffusion induced by red cell tumbling, it is possible that this mech- 
anism is also responsible for the observed aggregate geometry. An aggre- 
gate's growth is enhanced at its peak and reduced in the valleys between 
(3) 

aggregates , due to the dependence of platelet "diffusivity" on the local 
velocity gradient. Thus, the aggregates steepen as they grow, until the 
velocity gradient and diffusivity vanish in the valleys at a critical value of 
the ratio of aggregate size to separation distance. For longer times, an 
aggregate surface area with positive platelet flux is constantly reduced, and 
its ultimate size is soon attained. For a given aggregate half-wavelength X, 
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there is a characteristic time when the aggregate "sharpens" and achieves 
its maximum height. To estimate this time, consider the approximate 
solution^ 3 ) for the velocity gradient distribution on the wavy wall represen- 
ting an array of aggregates shown in Fig. 24. 


( — ) 

1 9z J w 



(34) 


When ^ 1/2 tt, the velocity gradient vanishes in the valleys and is 

enhanced by a factor of two at the peaks. As the aggregate grows, the 
growth rate increases at the peaks, since the diffusivity is given by 

D w = Cm2 / Sec 

and the platelet concentration gradient, for constant particle flux between 
streamlines, is given by 


(~) ~( J ( illL ) ~Jr 
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Therefore, a first-order approximation- of the platelet flux to the wall is 


F = D (-?-^-) 
w w 9 z 
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The aggregate growth rate is then approximated by: 
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Fig. 24 Schematic of platelet aggregate array growth model. 
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_ 1 2 3 

where Vp is the volume of a platelet, ~ 4 x 10 cm . Since the above 
equation is valid to first-order for t > 0, the following solution is found: 


z = Xt + e exp [3*0 sin ~ 
w o c L J X 


xT = 


with e o the amplitude at t = o. 

The first term in the equation above represents uniform growth of the 
aggregates, while the second term describes the unstable growth of the 
peaks. It is the second term that limits the growth of an aggregate, 
and this occurs on the time scale given by 


t ~ 1/3* 

The characteristic growth time is therefore 

“l ~ 1 

T G = X 3* BV (f^- ) (|£- ) 

G L P 9 Z co 

' 

Thus, the time to grow a mature aggregate of wavelength 2X is Tq. 
Conversely, the maximum aggregate growth in a time Tq is to be expected 
for aggregates separated by a distance of 2X, defined above. Denoting 
this distance as 2Xq, it follows that for X < Xq growth stops at t < Tq, 
while for X > Xq the growth rate is less than that for Xq. Thus, with 
increasing time, the initial large number of small closely spaced aggregates 
should be replaced by fewer large aggregates, with greater separation. 

While the above conclusion is in qualitative agreement with 
experimental observations, it remains to be shown that the calculated 
platelet flux, when used in Eq. (39), is consistent with the growth of 
a gg r ©g a te height and interaggregate distance with time. For a sufficiently 
larger shear rate, e. g. , outside the white cell circle, the Newtonian flow 
field solution^ 3 ) may be used. The unperturbed platelet flux then becomes: 
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o 3 _ 2 

Taking c^ = 4. 6 x 10 platelets/cm , the flux at A = 9 sec , = 0. 9 mm 

is 8 x 104 platelets/cm^ - sec and the expected interaggregate distance is 


2 Xq « . 05 r Q(microns) (41) 

To compare this calculation with experimental data, the 
distribution of interaggregate distances (taking all adjacent aggregate 
pairs) was measured at three times in a thrombus formed at the above 
conditions. This data is presented in Fig. 21, with a "best fit" curve 
drawn through each set. The interval of 2 \ with the maximum number 
of aggregates is taken as the characteristic interaggregate distance, and 
is compared. with 2 Xq in Fig. 25. The data are seen to be in remarkable 
agreement with the theoretical result, particularly in view of the many 
approximations involved in the analysis. It should also be noted that 
there appears to be, as expected, a slight dependence of 2X on radial and 
angular location in the thrombus. The latter effect is demonstrated by the 
double-peaked distribution in Fig. 21 which reflects the "age" of 
aggregates along the azimuthally growing edges of the thrombus. 

Two additional observations strengthen the model of aggregate 
growth described. First, the other possible growth mechanism, that of 
direct interception of platelets by aggregates, yields a growth rate that 
is an order-of-magnitude smaller than the diffusion rate. Second, the 
ability of platelet aggregates to trap red cells in their valleys depends on 
development of low velocity, low shear regions that are approximately 
an order-of-magnitude larger than the red cell diameter (8j Since red 
cell trapping (and reverse flow in some cases) is observed, and at times 
when both calculated and measured to be in the ~50 jj. range, there is little 
doubt that the aggregate and its surrounding flow field behave in the 
manner described. 
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Fig. 25 Interaggregate distance as a function of time. 
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In summary, it appears that platelet aggregates grow by an unstable 
diffusion/surface shape interaction, interaggregate distances scale with 
the local platelet flux rate and time, and the growth pattern is responsible 
for red cell trapping. Further investigations should consider the effects 
of three-dimensional aggregate arrays, more general growth patterns, 
and specific platelet aggregation mechanisms. 
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VII. CONCLUDING DISCUSSION 


The experimental data we have shown for the shape, size and 
growth rate of a thrombus on a foreign surface clearly show the effects 
of the flow of the blood on thrombus growth. Our understanding of the 
physical processes affecting the evolution of a thrombus, though hindered 
by the lack of knowledge of the detailed chemistry, has been increased 
by the application of simple fluid mechanical and mass transfer models. 

The basic concept underlying the growth model is the existence of a 
balance between the rate at which an activating species diffuses and is 
convected from the thrombus and the rate at which new platelets can 
diffuse from the flowing -blood to the surface. The flow dependence of this 
platelet diffusion rate is emphasized as the growth rate limiting factor, 
and the diffusion and convection of an activating species is shown to control 
the shape. 

The asymptotic character of the thrombus growth is treated and 
succeeds in demonstrating the flow dependence of the size and spacing of 
aggregates, as well as their role in trapping the larger red cells. Though 
the effect of the type of- surface material is not considered, we note the 
role of a surface inhomogeneity in triggering the initial platelet aggregation 
leading to the thrombus. The platelet aggregation is shown to have an 
unstable "switch-on" nature which occurs when the critical concentration 
of activating species is exceeded. It seems that a smooth artificial 
surface sufficiently free of surface imperfections may be quite non- 
thrombogenic under some flow conditions even if it does not prevent clot 
formation of stagnant blood. To settle questions of this type, further 
experimentation and mathematical modeling is certainly needed. 
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